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Summary 



Infection of bovine leucocytes by Theileria annulata 
results in establishment of transformed, infected 
cells. Infection of the host cell is known to promote 
constitutive activation of pro-inflammatory tran- 
scription factors that have the potential to be 
beneficial or detrimental. In this study we have 
compared the effect of LPS activation on uninfected 
bovine leucocytes (BL20 cells) and their Thellerla- 
infected counterpart (TBL20). Gene expression 
profiles representing activated uninfected BL20 
relative to TBL20 cells were also compared. The 
results show that while prolonged stimulation with 
LPS induces cell death and activation of NF-kB in 
BL20 cells, the viability of T'/re/Ver/a-infected cells 
was unaffected. Analysis of gene expression net- 
works provided evidence that the parasite estab- 
lishes tight control over pathways associated with 
cellular activation by modulating reception of 
extrinsic stimuli and by significantly altering the 
expression outcome of genes targeted by infection- 
activated transcription factors. Pathway analysis 
of the data set identified novel candidate genes 
involved in manipulation of cellular functions asso- 
ciated with the infected transformed cell. The data 
indicate that the T. annulata parasite can irrevers- 
ibly reconfigure host cell gene expression networks 
associated with development of inflammatory 
disease and cancer to generate an outcome that 



Received 25 January, 2012; revised 29 March, 2012; accepted 
19 April, 2012. *For correspondence. E-mail brian. shiels® 
glasgow.ac.uk; Tel. (+44)(0) 141 330 5756; Fax (+44)(0) 141 330 
5603. 

tpresent address: Wellcome Trust Centre for Gene Regulation 
and Expression, College of Life Sciences, University of Dundee, 
Dundee, UK. 

©2012 Blackwell Publishing Ltd 



is beneficial to survival and propagation of the 
infected leucocyte. 

Introduction 

The protozoan parasites Theileria annulata and T. parva 
cause economically innportant lynnphoproliferative disor- 
ders of cattle in developing countries of the Old World 
(Irvin and Morrison, 1987). In the laboratory they provide 
a unique model to investigate how a simple eukaryotic cell 
can direct gene expression and induce proliferation of 
another, more complex, eukaryotic host cell. Following 
invasion of bovine leucocytes by tick-inoculated sporo- 
zoites, rhoptries and microspheres discharge and the 
enveloping host cell membrane is disrupted. The interme- 
diate trophozoite stage is rapidly surrounded by an array 
of host microtubules and within 48 h differentiates into the 
macroschizont stage (Shaw and Tilney, 1 992). This differ- 
entiation event is accompanied by parasite-dependent 
transformation of the host cell (reviewed in Dobbelaere 
and Heussler, 1999), which in vivo manifests in exponen- 
tial proliferation of infected leucocytes over 4-5 days, 
followed by dissemination to a wide range of tissues. In 
acute cases of disease, disorganization and destruction of 
the lymphoid system occurs and death is often caused by 
massive pulmonary oedema associated with migration of 
infected cells to the lungs (Irvin and Morrison, 1987). 

Transformation of T. annulata- and T. pa/va-infected 
leucocytes is known to be dependent on constitutive acti- 
vation of bovine transcription factors and associated 
signal transduction pathways (reviewed in Dobbelaere 
and Heussler, 1999; Shiels etal., 2006). Activation of 
AP-1, ATF-2, cMyc, NF-kB and STAT3 have all been 
documented and are likely to be responsible for many of 
the phenotypic alterations displayed by the infected cell. 
However, in addition to being beneficial, constitutive tran- 
scription factor activation has the potential to be detrimen- 
tal to survival of the infected cell or host. For example, 
constitutive activation of c-Jun (a dominant component of 
AP-1 transcription factor), is thought to play an important 
role in metastasis of TtoVer/a-infected cells (Adamson 
etal., 2000), while activation of NF-kB confers protection 
against apoptosis (Heussler etal., 1999). However, both 
these factors can be pro-apoptotic (Dunn etal., 2002; 
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Dutta et al., 2006) by upregulating genes encoding death- 
inducing signals, including pro-inflammatory cytokines. It 
has been predicted that NF-KB-dependent regulation of 
inflammatory genes contributes to the pathogenesis of 
theileriosis (Heussler etal., 2002). 

A recent study demonstrated that only a proportion of 
leucocytes infected with T. parva sporozoites successfully 
progress to proliferation, while the rest undergo apoptosis 
(Rocchi etal., 2006). Thus, events that are detrimental 
to the infected cell can occur following sporozoite inva- 
sion and transformation is not guaranteed. In addition, 
although basal expression of inflammatory/cellular 
defence genes such as TNFa, iNOS and ISG15 can be 
detected (Sager etal., 1997; Oura etal., 2006), Theileria- 
infected cells are refractory to LPS-mediated stimulation 
of high-level expression, despite constitutive activation of 
transcription factors (e.g. AP-1 and NF-kB) associated 
with this event. Furthermore, parasite-mediated repres- 
sion of highly abundant mRNA levels has been indicated 
for the ISGylation pathway genes, UBP43 and ISG15 
(Oura etal., 2006) whose expression, like TNFa and 
INOS, can be induced by activated NF-kB (Li et al., 2001 ; 
Liu etal., 2011). Taken together, these studies predict that 
Thelleria parasites must direct gene expression networks 
associated with activation of host cell transcription factors 
to promote an outcome that favours establishment of the 
transformed cell. 

In this study we have conducted a comprehensive analy- 
sis of expression profiles associated with LPS stimulation 
of the uninfected bovine lymphosarcoma cell line BL20 
(Olobo and Black, 1989) and the 7^ annu/afa-infected 
counterpart, TBL20 (Shiels etal., 1986). These lines were 
chosen as they provide uniform cellular populations to 
allow robust detection of differences in gene expression 
associated with altered cellular phenotypes. Moreover, 
TBL cells have been used to investigate infection- 
associated phenotypes (Shiels etal., 1986; Baylis etal., 
1 995) and are known to be dependent on a viable parasite 
for survival, as they undergo apoptosis when the macro- 
schizont is killed by BW720c (buparvaquone) (Dobbelaere 
and Heussler, 1999). The results of our analysis indicate 
that while the modulated gene expression profiles of LPS- 
stimulated, uninfected BL20 and infected TBL20 cells sig- 
nificantly overlap, a major reorganization of activated gene 
expression is orchestrated by the parasite-infected cell to 
promote its establishment and alter receptiveness to 
extrinsic activation signals. 

Results and discussion 

Differential growth and activation response of BL20 and 
TBL20 cells to stimulation by LPS 

From the results of previous studies we postulated that 
uninfected cells could display a detrimental response to 

©2012 Blackwell Publishing Ltd, Cellular Microbiology, 14, 1434-1454 



Modulation of host gene expression by J. annulata 1 435 

prolonged stimulation with an inflammatory mediator, 
while The/Ver/a-infected cells would be refractory to such 
stimulation. To test this postulation, BL20 and TBL20 cul- 
tures were treated with LPS and the number of viable and 
dead cells estimated after 18 h. As shown by Fig. 1A, the 
number of viable BL20 cells following stimulation with 
LPS was significantly lower (2.3-fold) relative to the 
un-stimulated control. This loss of growth potential was 
accompanied by a 4-fold increase in the percentage of 
non-viable cells (Fig. IB) and a 2.9-fold increase in 
caspase 3/7 activity (Fig. 1C). Infected TBL20 cells, in 
contrast, showed no significant difference between LPS- 
stimulated and control cultures. These results suggest 
that TBL20 cells are resistant to a detrimental effect of 
stimulation with LPS (18 h). 

In addition to growth analysis, infected and uninfected 
cells were assayed for NF-kB reporter activity. NF-kB 
reporter activity was significantly elevated in uninfected 
BL20 cells following LPS stimulation and, while basal 
reporter activity was higher in TBL20, no significant eleva- 
tion in constitutive NF-kB activity following LPS treatment 
was detected (Fig. 1C). Confirmation of NF-kB activation 
by LPS in BL20 was obtained by detection of the p65 
subunit of NF-kB in the nucleus and IKKy signalosome 
complexes in the cytoplasm (Fig. SI). The finding that 
NF-kB activation occurs in LPS-stimulated BL20 and is 
constitutively activated in infected TBL20 suggests that for 
this transcription factor, at least, TBL20 cells may avoid 
a potentially detrimental response by modulating the 
outcome of the activation event. 

Activation of BL20 cells by LPS gives rise to gene 
expression changes that show significant similarity to 
changes associated with infected TBL20 cells 

To compare bovine gene expression changes associated 
with activation of BL20 cells with LPS relative to changes 
associated with infected TBL20 cells, microarray analysis 
was performed. The oligonucleotide microarray platform 
represented all confirmed and putative bovine mRNA fea- 
tures available (31 481) and provided an increase in gene 
number relative to the BoMP array utilized by Jensen 
etal. (2008). A subset of 19 777 genes was used for 
further analysis, as this represents sequences on the 
array that were found to be present in the current bovine 
genome assembly or RefSeq database. The array was 
hybridized with cDNA representing RNA from cultures of 
BL20 and TBL20 cells before and at 4 and 18 h post LPS 
stimulation. Differential gene expression on the normal- 
ized data set was assessed by pair-wise comparison 
using Rank Product Analysis (Breitling etal., 2004) and 
differences were considered to be significant at a false 
discovery rate (FDR) of less than 5% and a fold change 
greater than 2. 
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Fig. 1. A. Viability of BL20 and TBL20 cells following LPS stimulation. Cell viability of BL20 and TBL20 cell cultures was assessed 

using Trypan Blue exclusion in the presence of either LPS (1 |ig ml"') or no drug DMSO (1 \i\ ml"') vehicle control (ND) at 18 h 

post treatment. Counts are expressed as group mean ± SD, n = 3. The number of viable cells in BL20 LPS-treated cultures was significantly 

lower than BL20 no drug (ND) and TBL20 cultures (f-test, P< 0.01); TBL20 compared with TBL20-LPS cultures showed no significant 

difference. 

B. Percentage dead cells in BL20 and TBL20 cultures following LPS stimulation. Cell viability of BL20 and TBL20 cell cultures was assessed 
using Trypan Blue exclusion in the presence of either LPS (1 |ig ml"') or no drug (ND) at 18 h post treatment. The percentage dead cells in 
each culture condition are expressed as group mean ± SD, n = 3. The percentage dead cells in BL20 LPS-treated cultures was significantly 
greater than BL20 no drug (ND) and TBL20 cultures (f-test, P< 0.01); TBL20 compared with TBL20-LPS cultures showed no significant 
difference. 

C. Caspase 3/7 activity in BL20 and TBL20 cells following LPS stimulation. Caspase 3/7 activity was determined after incubation of cells in 
the presence of either LPS (1 \ig ml"') or the no drug DMSO (1 |il ml"') vehicle control (ND) for 18 h. Results are expressed as relative 
fluorescence units (mean ± SD, n = 3). Caspase activity was significantly greater in BL20-LPS compared with BL20 no drug (f-test, P < 0.01 ); 
TBL20 compared with TBL20-LPS cultures showed no significant difference. 

D. NF-KB-dependent reporter activity in BL20 and TBL20 cells following LPS stimulation. Levels of NF-kB activation in Theileria annulata- 
infected TBL20 cells and uninfected BL20 cells following stimulation with LPS are shown relative to untreated controls. Cells were transiently 
co-transfected with an NF-kB promoter/Firefly luciferase reporter construct and a control Renilla luciferase construct. Twenty-four hours after 
transfection, cells were stimulated either with LPS (1 |ig ml"') or no drug DMSO (1 |il ml"') vehicle control (ND) for 4 h before extraction and 
assay for luciferase activity. NF-kB promoter activity is expressed as the ratio of NF-kB Firefly luciferase/Renilla activity (mean ± SEM, n = 4). 
NF-kB activity was significantly higher in BL20-LPS compared with BL20 no drug (f-test, P< 0.005); TBL20 compared TBL20-LPS cultures 
showed no significant difference. 



The results showed that TBL20 cells were refractory to 
major LPS-induced global expression changes, as only a 
small subset of genes (214 genes) showed significantly 
altered expression following stimulation of TBL20 with 



LPS, whereas a total of 2180 genes were identified as 
altered by LPS stimulation of BL20 (BL20-LPS data set). 
Thus, as implied from previous studies, T. annulata- 
infected leucocytes are associated with major suppression 
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Fig. 2. Identification of a significant overlap between TBL20(IA) 
and BL20-LPS data sets. Analysis of the lists of differentially 
expressed genes (FDR < 5%, FC > 2) generated from each 
experimental condition led to the identification of a set of 794 
genes common to data sets representing parasite-infected cells 
and LPS treatment of uninfected cells. The overlap between two 
data sets was found to be highly significant (P< 0.001). 

of lipopolysaccharide (LPS) inducible gene expression 
(Sager etal., 1997; Oura etal., 2006; Jensen etal., 2008). 

A large number of genes (1811) showed changes in 
expression that were associated with infection of BL20 
cells by T. annulata [the infection-associated data set, 
TBL20(IA)]. Comparison of the BL20-LPS data set and 
the TBL20(IA) data set showed the common overlapping 
gene list (794 genes) was significantly larger than would 
be expected to occur by chance (Fig. 2). This result dem- 
onstrates that even though it is largely unresponsive to 
LPS, the gene expression profile of the parasite-infected 
cell is related to an activated state induced by a classical 
inflammatory mediator. Pathway analysis of the BL20- 
LPS data set confirmed that enrichment for inflammatory 
response genes was of greatest significance (data not 
shown). As infection of the host leucocyte by T. annulata 
results in constitutive activation of transcription factors 
such as NF-kB (Oura etal., 2006) and AP1 (Baylis etal., 
1995) that are important regulators of the inflammatory 
response, it was of interest to investigate whether the 
outcome of their activation is altered in the parasite- 
infected cell and results in a prediction of phenotype that 
is favourable to the establishment of infection. 

Gene expression changes in response to activation 
of BL20 cells by LPS are modulated in 
T. annulata-/nfecfec/ TBL20 cells 

To investigate whether gene expression changes associ- 
ated with cellular activation of BL20 cells are modulated in 

©2012 Blackwell Publishing Ltd, Cellular Microbiology, 14, 1434-1454 
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parasite-infected cells, gene expression in TBL20 cells 
was contrasted with that of uninfected, activated cells. To 
achieve this, the set of differentially expressed genes in 
the BL20/BL20-LPS data set was compared against the 
expression values from the entire TBL20 data set. For an 
LPS stimulated gene to be considered as a candidate for 
modulation by infected cells, the expression level in the 
infected cell had to show a difference of > 2-fold com- 
pared with that of the LPS stimulated BL20 cell (BL20- 
LPS: TBL20 ratio > 2), or show a contrasting Rank 
Products result. This resulted in the generation of an LPS 
response, infected cell-modulated gene list of 1214 genes 
(hereafter known as the infected cell-modulated gene list). 
Of these, 1 98 were more highly expressed in TBL20 cells 
than in BL20 LPS-treated cells, while a considerably 
larger number, 1016 genes, were repressed in TBL20 
cells. As expected (see above), the vast majority (> 95%) 
of the genes in the list of 1214 showed no significant 
change in their expression level in TBL20 cells stimulated 
with LPS (FDR < 5%). 

Further analysis of the data set was conducted by 
incorporating data obtained from an independent study 
on the response of host cell gene expression following 
treatment of TBL20 with BW720c for 48 h (J. Kinnaird 
etal., in preparation). This enabled identification of 
expression changes that are more likely to be under 
direct control of a viable parasite, i.e. if an infection- 
associated modulation was reversed upon BW720c 
treatment (see Tables S1-S8). In addition, analysis of 
expression values for the set of 1214 genes across all 
six cellular conditions partitioned the gene list into four 
pairs of reciprocal profiles with evidence of modulated 
expression in TBL20 relative to BL20-LPS (Figs 3 and 4 
for selected genes and Tables S1-S8 for full gene lists 
of each profile). To provide insight into potential biologi- 
cal relevance of the infected cell-associated modulation, 
the data set for each profile was subjected to pathway 
analysis to investigate enrichment for predicted molecu- 
lar and cellular functions. Functional annotation derived 
from the bovine RefSeq database was also utilized. For 
validation of the microarray results, between three and 
nine genes were chosen from each profile for semi- 
quantitative reverse transcription PGR (RT-PCR). These 
were selected based, primarily, on annotation that could 
provide possible insight into the infected cell modulation 
event. Two validated genes of interest were then 
selected for further validation by qRT-PCR. 

Profiles 1 and 2 represent an enhancement of gene 
expression in infected cells compared with the modulation 
detected for LPS stimulated BL20 (elevation of an LPS 
upregulated gene, profile 1 ; stronger repression of an LPS 
downregulated gene, profile 2). These patterns of modu- 
lation could reflect a quantitative difference in activation 
status in the infected cells. The level of NF-kB reporter 
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Profile 



Gene 
Symbol 



Annotation 



Modul- Example RT PCR results 

ation BL20 lps tbl2olps 

18h 4h Oh Oh 4h 18h 



'Super-up In 
Infected cell' 



BL20 BL20 TBU0TBL20 
LPS LPS 



SLAMF7 SLAM family 7; 19A protein 

LRRTM4 Leucin rich repeat transmembrane neuronal 4 

MMP13 Matrix metallopeptidase 13 

XCL1 Chemokine (C motif) ligand 1 

TP73 Tumor protein P73 

IL2RB Interleukin 2 receptor, beta 

ICAM1 Intercellular adhesion molecule 1 

SELL Selectin L 

MEIS3 Meis homeobox 3 

SLC6A6 Solute carrier family 6 

CTTN Cortactin 



+148 



(n= 32) 



+68 






MMP13 


+42 


MEIS3 


+20 




+ 11 


CTTN 


+ 12 




+ 10 


IL2RS 


+9 




+8 


SELL 


+4 






ICAIVI1 




SLC6A6 


-88 


MYL4 


-85 




CD36 


-71 




-44 


CALD1 


-26 




-17 


CCL3 


-14 




-11 


IVIFAP2 


-9 
-8 


LIVI02 


-6 


PGM5 


-4 


IL18 




TLR4 




'Super-down in 
Infected cell' 



BL20 BL20 TBL20 TBL20 
LPS LPS 



(n= 180) 



LOC404176 T-cell receptor delta chain 

IGF1R Insulin-like growth factor 1 receptor 

F13A1 Coagulation factor XIII A chain 

TLR4 Toll-like receptor 4 precursor 

CCL3 Chemokine (C-C motif) ligand 3 

MYL4 Myosin, light chain 4 like protein 

PGM5 Phosphoglucomutase 5 

LM02 LIM domain only 2; rhombotin-like 1 

CALD1 Caldesmon 1 

MFAP2 Microfibrillar-associated protein 2 

CD36 CD36 molecule, thrombospondin receptor 

IL18 Interleukin 18 




3 

'Down In Infected 
cell' 



BL20 

LPS 



BL20 TBL20TBL20 
LPS 



(n = 154) 



MPEG1 


Macrophage expressed gene 1 


-449 


CLU 


CLU 


Clusterin 


-281 




CD200R 


CD200 receptor 1-like 


-133 


IL15 


BOLA-DYA 


MHC, class II, DY alpha 


-132 




EPHA7 


EPH receptor A7 tyrosine kinase 


-120 


STAP1 


IL15 


Interleukin 15 


-107 


IGJ 


BTN3A2 


Butyrophilin, subfamily 3, member A3 


-83 


IGJ 


Immunoglobulin J chain 


-58 


IVIPEG1 


STAP1 


Signal transducing adaptor family member 1 


-44 




BIRC3 


Baculoviral lAP repeat-containing 3 


-26 


BoLA-DYA 








BIRC3 




4 




AHNAK 


AHNAK nucleoprotein 


+235 


ANGPT4 


'Up in Infected cell' 


LOC514557 


OAF homolog (Drosophila) 


+134 




c 




CD99 


CD99 molecule; HBA71 


+113 


NDRG1 


o ^, 




LOC5 12831 


Heterogeneous ribonucleoproteins A2/B1 


+106 




tj ^ 




THEM4 


Thioesterase superfamily member 4 


+ 104 


PAK1 






SPARC 
NDRG1 


Secreted protein, acidic, cysteine-rich 
N-myc downstream regulated gene 1 


+53 
+43 




BL20 BL20 TBL20TBL20 


RNASEL 


LPS 


LPS 


ANGPT4 


Angiopoietin 4 


+36 








RNASEL 


Similar to 2-5A-dependent RNase 


+31 


PIK3R6 






PAK1 


p21 protein (Cdc42/Rac)-activated kinase 1 


+22 








BCL9 


B-cell CLL/lymphoma 9 


+21 




(n = 55) 




PIK3R6 


Phosphoinositide-3kinase, regulatory subunit 6 


+8 






RT PCR control 



ACTB 



Actin, beta 



ACTS 



^^^^^ 



Fig. 3. Expression profiles of BL20-LPS-response genes showing a modulated response in rte//eria-infected cells. The different gene 
expression response profiles for TBL20 and BL20 cells obtained before and following LPS treatment are shown diagrammatically to the left, 
together with the total number of genes in that profile. Denoted genes (middle) represent top-ranking/selected (based on predicted function) 
differentially expressed genes within the profile. Semi-quantitative RT-PGR validation for a panel of the selected genes is shown on the right. 
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Profile 



Gene 
Symbol 



Annotation 



Modul- 
ation 



Example RT PCR results 

BL20 LPS TBL20LPS 

18h 4h Oh Oh 4h 16h 



IE 

9 


OAS2 


2'-5'-oligoadenylate synthetase 2 


-32 


CD69 


'Up less in Infected 


CD83 


CD83 molecule 


-28 




cell* 


CXCL10 


Chemokine (C-X-C motif) ligand 10 


-24 


ISG16 


c 


RELB 


v-rel reticuloendotheliosis viral oncogene B 


-11 


CXCL10 


■~ *\ 


TNFAIP3 


Hypothetical 


-9 


'£ ^^--^ 

Q. 


TGFBR3 


Transforming growth factor, beta receptor III 


-8 


MX2 




CX3CR1 


Chemol<ine (C-X3-C motif) receptor 1 


-8 








BL20 BL20 TBL20TBL20 


TIFA 


TRAF-lnteracting protein, forl<head-ass. 


-8 


CX3CR1 


LPS LPS 


CX3CR1 


Chemokine (C-X3-C motif) receptor 1 


-8 






TIFA 


Hypothetical 


-7 






CD69 


CD69 molecule 


-7 




(n = 38) 


ISG15 


ISG15 ubiquitin-like modifier 


-5 






'Down less In 
Infected cell' 



BL20 BL20 TBL20TBL20 
LPS LPS 



(n = 13) 



LOC528166 


Uncharacterlsed protein 


+5 


LOC781353 


Derailed-like RYK 


+4 


CYP3A4 


Cytochrome P450, 3A, 4 


+4 


STXBP6 


Syntaxin binding protein 6 


+4 


RYK 


RYK receptor-like tyrosine kinase 


+4 


CYP4F3 


Cytochrome P450, 4F, 3 


+3 


LOC528166 


Uncharacterlsed protein 


+3 


PGLYRP1 


Peptidoglycan recognition protein 1 


+3 


TRAF3IP3 


TRAF3 interacting protein 3 


+3 


LOC537817 


TBG1 domain family, member 9 


+3 


GALNT4 


N-acetylgalactosaminyltransferase 4 


+3 


AQP3 


Aquaporin 3 


+2 




7 


VCAM1 


Vascular cell adhesion molecule 1 


-185 


VCAM1 


'NC In Infected cell 


TFEC 


Transcription factor EC 


-47 




though up in LPS' 


LOC511753 


Olfactory receptor, family 5, subfamily 1 


-40 


CCL4 




LOC784595 


Sim. olfactory receptor Olfrl 174 


-27 




i \ 


ADAMTSL5 


ADAMTS-like protein 5 


-27 


RUNX2 


a > • • 


CCL4 


Chemokine (C-C motif) ligand 4 


-17 


IFIT3 


UJ 


IFIT3 


IFN-induced prot. withtetratricopeptide rpts 


-15 


BL20 BL20 TBL20 TBL20 




3 




DTX4 


LPS LPS 


DTX4 


DTX4 protein 


-11 




CFLAR 


CASP8 and FADD-like apoptosis regulator 


-9 


CFLAR 




RUNX2 


Runt-related transcription factor 2-like 


-5 




(n = 650) 



'NC In Infected cell 
though down In LPS' 



BL20 BL20 TBL20TBL20 
LPS LPS 



(n = 98) 



PLIN3 


Perilipin 3 


+17 


BCOR 


Sim. BCL-6 co-repressor 


+15 


SERPINE2 


Serpin peptidase inhibitor, E2 


+15 


FGD4 


FYVE, RhoGEF and PH domain containing 1 


+15 


CAMKV 


CaM kinase-like vesicle-associated 


+11 


ARMCX6 


Armadillo repeat containing, X-linked 6 


+10 


LOC616063 


Sim. myeloid-assoc. differentiation marker 


+10 


NID2 


nidogen 2 (osteonidogen) 


+9 


SLC7A11 


Solute carrier family 7 


+9 


LOC537873 


Uncharacterlsed protein 


+9 


LOC540836 


Uncharacterlsed protein 


+7 




RT PCR control 



ACTB 



Actin, beta 



ACTR 



LUJJJJ 



Fig. 4. Expression profiles of BL20-LPS-response genes showing a modulated response in rrte/teria-infected cells. The different gene 
expression response profiles for TBL20 and BL20 cells obtained before and following LPS treatment are shown diagrammatically to the left, 
together with the total number of genes in that profile. Denoted genes (middle) represent top-ranking/selected (based on predicted function) 
differentially expressed genes within the profile. Semi-quantitative RT-PGR validation for a panel of the selected genes is shown on the right. 
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activity, for example, was found to be significantly higher 
in TBL20 compared with BL20-LPS. 

Profile 1 comprised 32 genes and for many of these 
genes the amplitude of enhanced, expression associated 
with infected cell was considerable (Fig. 3 and Table S1). 
Thus, the genes encoding chemokine (C motif) ligand 1, 
MMP13 and an IL2RB-like protein were elevated 42-, 54- 
and 11 -fold respectively in TBL20 relative to the maximal 
response obtained for LPS stimulated BL20 cells. Of the 
14 BW720c-sensitive genes identified (Table SI), 10 
(71.4%) responded in an expected manner, i.e. expres- 
sion was significantly reduced by treatment over 48 h. 
Semi-quantitative RT-PCR was performed for eight 
genes, each of which validated the microarray profile 
(Fig. 3); the gene encoding the MEIS3 transcription factor, 
for example, showed a clear elevation in PGR product in 
LPS-stimulated uninfected BL20 cells but a considerably 
higher level in infected TBL20, with no obvious detectable 
alteration following LPS stimulation of TBL20. Further vali- 
dation of the profile was obtained by qRT-PCR for cortac- 
tin and MEIS3 (see Fig. S2, panels 1 A and 1 B) where the 
results agreed with the microarray data. qRT-PCR for 
ME/S3 showed a 5.98-fold difference (log2) in TBL20 rela- 
tive to BL20-LPS, compared with a difference of 3.1 -fold 
(log2) using the array data, while qRT-PCR for cortactin 
calculated the difference to be 1.31 -fold (log2) compared 
with 1.18-fold (log2) from the array data. The two func- 
tional categories that showed the most significant enrich- 
ment for profile 1 genes (Fig. S3) were 'control of cellular 
movement' and 'cell-to-cell signalling', and there was also 
an enrichment of genes associated with inflammation in 
this profile. In addition, a number of genes in profile 1 
have predicted functions that can be linked to the pheno- 
type of T/7e/7era-transformed leucocytes. Thus, elevated 
expression was detected for a matrix metallopeptidase 
gene (Baylis etal., 1995), the IL2 receptor B-like gene 
(Herrmann etal., 1989), genes encoding cellular adhe- 
sion molecules, ICAM1 (Jensen etal., 2008) and selectin 
L (SELL), and genes that regulate the cytoskeleton, cort- 
actin (CTTN) (Baumgartner, 2011). It can be postulated 
that genes captured in profile 1 are expressed at an 
elevated level in TBL20 because they are advantageous 
for establishment and or survival of the infected cell. 

Profile 2 is the reciprocal of profile 1 and represents 180 
genes whose expression is repressed in BL20-LPS but 
super-repressed in TBL20. Twenty-four of these genes 
showed major repression of expression associated with 
the infected cell (> 10-fold). Of the 83 genes sensitive to 
BW720C, 99.8% responded in the expected manner, indi- 
cating repression requires a viable parasite. Nine genes 
of this profile were validated by RT-PCR and agreed with 
the array profile (Fig. 3). qRT-PCR validation was per- 
formed for TLR4 and MYL4 and showed good correlation 
with the microarray results (Fig. S2, panels 5Aand 5B). A 



number of genes encoding transcription factors were 
present in the data set including LEF-1 (down 15-fold) and 
LM02 (down 11 -fold), both of which are strongly associ- 
ated with leukaemia (Curtis and McCormack, 2010; Guti- 
errez etal., 2010). Moreover, genes encoding membrane 
proteins/receptors were particularly prominent in the data 
set. Among 13 receptor genes that showed stronger 
downregulation in TBL20 compared with BL20-LPS were 
those encoding the IGFR1 receptor (down 85-fold in 
TBL20 relative to BL20-LPS), the Toll-like receptor 4 
(TLR4: down 44-fold), the G-protein-coupled receptor 
(GGR6: down 23-fold), S1PR1 (down 6-fold) and the 
thrombospondin receptor (down 6-fold). Thus, data from 
this profile strongly suggests that genes encoding mol- 
ecules located to the surface of the infected cell are 
repressed and could influence interaction of the infected 
cell with a range of ligands, inflammatory mediators and 
cells of the bovine immune system. It can be predicted 
from the super-repression profile that high-level expres- 
sion of the identified genes is most likely to be detrimental 
to the infected cell. 

TLR4 is the major cellular receptor that transduces the 
LPS activation signal (Chow etal., 1999) and significant 
reduction in expression (61 .4-fold between BL20 versus 
TBL20) could play a role in the unresponsiveness of 
777e/7er/a-infected cells to LPS identified in this and previ- 
ous studies. The gene encoding the LPS ligand GDI 4 was 
also expressed at a very low level in both uninfected BL20 
and TBL20 cells. It has been demonstrated, however, that 
a lack of GDI 4 expression in bovine cells does not pre- 
clude recognition of LPS stimulation, possibly via the 
presence of soluble GDI 4 present in fetal calf serum 
(Sauter etal., 2007). 

Super-repression of TLR4 was validated at the protein 
level: immunoblotting with a specific anti-TLR4 antibody 
clearly demonstrated lower levels of TLR4 protein in 
TBL20 extracts relative to uninfected BL20 extracts from 
control or LPS stimulated cells (see Fig. 5A). In addition, 
both RT-PCR and immunoblotting showed that killing the 
parasite with BW720c considerably elevated the level of 
TLR4 mRNA and protein (Fig. 5B and G). Thus, repres- 
sion of TLR4 by the infected cell is dependent on a viable 
parasite. The TLR4 signal transduction pathway plays a 
major role in activation of NF-kB following recognition of 
pathogen-associated molecular patterns (PAMPs). Sim- 
plistically, repression of this pathway may not be in the 
interest of the r;f7e//er/a-infected cell, since it has been 
demonstrated that constitutive activation of NF-kB is 
required for its survival (Heussler etal., 1999). Parasite- 
dependent constitutive activation of NF-kB, however, is 
considered to be independent of extrinsic signalling, as 
the parasite promotes aggregation of active IKK com- 
plexes at the macroschizont surface (Heussler etal., 
2002). Whether the repression of TLR4 gene expression 
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Fig. 5. A. Immunoblot analysis of TLR4 protein expression 
following stimulation of BL20 and TBL20 cells with LPS. Expression 
of TLR4 protein was detected by immunoblot analysis of total 
protein extracts derived from uninfected and The/ter/a-infected cell 
cultures stimulated with either LPS for 4 and 18 h or no drug 
DMSO (1 \i\ mM) vehicle control (ND). Detection of p-tubulin was 
used as a constitutive control to ensure equal loading of extract. 

B. RT-PCR analysis of TLR4 expression following treatment of 
TBL20 cells with BW720c. Expression of TLR4 at the mRNA level 
was assessed by semi-quantitative RT-PCR using RNA isolated 
from BL20 and TBL20 cells cultured with no drug (ND) or TBL20 
cells cultured with BW720c at 50 ng ml"' for 48 h. Primers specific 
for the p-actin gene (ACT P) were used as a control for constitutive 
expression. 

C. Immunoblot analysis of TLR4 protein expression following 
treatment of TBL20 cells with BW720c. Expression of TLR4 protein 
was assessed by immunoblot analysis of total protein extracts of 
BL20 and TBL20 cells cultured with no drug (ND) or TBL20 cells 
cultured with BW720c at 50 ng ml"' for 48 h. Detection of p-tubulin 
was used as a constitutive control to ensure equal loading of 
extract. 



provides an advantage to the infected cell is not clear but 
it could reduce the level of activation of NF-kB that occurs 
via extrinsic inflammatory stimuli. This would allow the 
level of NF-kB activation to be primarily set by the 
parasite-dependent IKK signalosome pathway. A block in 
TLR signalling has been reported for Toxoplasma gondii. 
It has been proposed that this allows prevention of pro- 
inflammatory responses triggered by T. gondli-demed 
TLR ligands, or gut flora, or reflects a general non- 
responsiveness to pro-inflammatory stimulation (reviewed 
by Leng etal., 2009). In addition, mammalian p38 MAPKs 
which are activated by LPS and play a major role in 
apoptosis (Guha and Mackman, 2001; Zhou etal., 2007) 
are known to have low or undetectable activity in 
T"he/7er/a-infected cells (Botteron and Dobbelaere, 1998). 



While modulation of expression associated with profiles 
1 and 2 may relate to the level of transcription factor 
activation, this is unlikely to be the case for profiles 3 and 
4. In these profiles TBL20 was associated with a gene 
expression change opposite to that observed for BL20- 
LPS, and it can be postulated that it is of benefit to the 
infected cell to strongly repress (profile 3) or elevate 
(profile 4) expression of these genes. Profile 3 represents 
154 genes elevated by LPS stimulation of uninfected cells 
that are significantly repressed in infected TBL20 (Fig. 3). 
Almost all of the 63 genes that were sensitive to BW720c 
treatment in this profile responded in a reversible manner 
(i.e. 98.4% showed elevated expression with BW720c 
treatment), indicating active repression by a viable para- 
site. Many genes (n = 84) provided evidence of a repres- 
sion corresponding to a fold change of greater than 
ten. Seven genes from profile 3 were validated by 
semi-quantitative RT-PCR (Fig. 3). The gene encoding 
MPEG1, for example, showed an elevation in PGR 
product in LPS-stimulated uninfected BL20 cells while the 
level was clearly lower in infected TBL20 relative to BL20 
cells. qRT-PCR performed for the BOLA-DYA and BIRC3 
genes (Fig. S2, panels 2A and 2B) validated the microar- 
ray profile. Profile 3 showed enrichment for genes pre- 
dicted to encode proteins involved in cell-to-cell signalling, 
cell morphology, cellular development and cellular growth 
and proliferation. A significant number of genes were 
associated with cancer and apoptosis. Genes in profile 3 
included those encoding macrophage expressed MPEG1 , 
the inflammatory cytokine I LI 5, and the major histocom- 
patibility complex class II DY alpha. Repression of these 
genes may be advantageous to the parasite-infected cell 
by altering the differentiation status of the host cell or 
manipulating the immune response. 

Profile 4 is the reciprocal of profile 3 and represents 55 
genes that were repressed by LPS treatment of unin- 
fected BL20 cells but were elevated in infected TBL20 
cells. The level of modulation associated with this profile 
was considerable. For example, the expression of genes 
encoding AHNAK, CMRF-35 and SLA were amplified 235- 
fold, 99-fold and 45-fold respectively. A further 33 genes 
were elevated more than 10-fold. Interestingly, a large 
proportion (61.5%) of the 13 BW720c-sensitive genes 
responded in an unexpected manner, i.e. they were 
upregulated by BW720c. Microarray expression data for 
profile 4 was validated by RT-PCR using five genes. qRT- 
PCR validation was performed for ANGPT4 and NDRG1 
and showed good correlation with the microarray (see 
Fig. S2, panels 6A and 6B). Of the 51 genes that could 
be mapped onto the pathway analysis framework, 21 
showed an association with cancer. The profile 4 data set 
showed enrichment for genes encoding proteins that 
function in 'cell movement', 'cell death' and 'cell growth 
and proliferation'. Gene ontology indicated that a number 
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of genes (ANGPT4, NDRG1, PAK1, PI3KR6, SLA, 
SPARC) operate in signal transduction events. Overall, 
this profile suggests tliat the infected cell is strongly asso- 
ciated with upregulation of genes linl<ed to neoplasia. 

Profiles 5 and 6 both provide evidence that infected 
TBL20 cells attenuate the response of uninfected BL20 
cells to LPS. Thus, profile 5 represents 38 genes with 
expression levels elevated in 7/?e//er/a-infected cells, but 
to a lesser extent than in LPS stimulated BL20 cells. An 
example of this profile is the gene for interferon-induced 
protein 15, ISG15, elevated 121 -fold in BL20-LPS but only 
26-fold by infection TBL20 (relative to BL20) and no 
change in TBL20-LPS. These l<inetics mirror reduced 
expression associated with 7he//er/a-infected cells previ- 
ously observed by Northern blotting (Oura etal., 2006). 
Five genes of the profile were validated by semi- 
quantitative RT-PCR (Fig. 4). qRT-PCR performed for the 
CD69 and ISG15 genes confirmed the microarray data 
(Fig. S2, panels 3A and 38). Functional analysis showed 
enrichment for genes encoding proteins involved in 
protein modification, the antimicrobial response, the 
inflammatory response, cellular defence and the innate 
immune response (e.g. ISG15). Genes encoding tran- 
scription factors activated by Theileria infection such as 
NF-kB2, NF-kBIZ and Jur) were also identified. Altered 
expression of genes in profile 5 could indicate a limited 
upregulation of a beneficial gene that does not occur to 
the same degree as in BL20-LPS (e.g. the Jun proto- 
oncogene). However, it may be more lil<ely that high-level 
expression of cellular defence genes, such as ISG15 
(Oura et al., 2006) that occur via activation of the cell, are 
repressed by infection as they are potentially detrimental. 
The enrichment in genes associated with the interferon/ 
innate immune/cellular defence response, e.g. 0AS2, 
ISG15, MX2, IRF5 and //./(Moller etal., 1996; Jin etal., 
1999; Urosevic, 2003; Pitha-Rowe and Pitha, 2007; Bin 
etal., 2011), supports this possibility. Furthermore, 72% of 
the 18 BW720c-sensitive genes in profile 5 showed 
elevated expression upon drug treatment, indicating 
active repression by the viable parasite. 

Profile 6 was the reciprocal of profile 5 and repre- 
sented genes where expression was repressed with LPS 
treatment of BL20 but repressed to a lesser extent in 
parasite infected TBL20 cells. This profile represented 
the smallest number of genes (n= 13), and four out of 
the five BW720c-sensitive genes were elevated upon 
treatment. The array profile was validated for three 
genes by RT-PCR and by qRT-PCR for CYP4F3 and 
RYK (Fig. S2, panels 7A and 7B). Pathway analysis of 
the data set for this profile showed enrichment for genes 
predicted to function in metabolism or small molecule 
biochemistry. 

Profiles 7 and 8 represent no change in infected TBL20 
cells compared with an altered response identified for 



BL20-LPS (up or down). Profile 7 represents 650 genes 
elevated by LPS stimulation of BL20 that showed no 
change on comparison of BL20 with TBL20 and was the 
largest obtained by the study. Only a small number of 
these genes {n = 6) were sensitive to BW720c, supporting 
a mechanism of non-responsiveness, rather than repres- 
sion of an activation event associated with the infected 
cell. In some cases, however, postulation of a lack of 
response mechanism may be over-simplistic. VCAM1, for 
example, is a classical NF-KB-dependent gene (Luo etal., 
2010) and it is conceivable that infection of the host cell is 
associated with a more direct repression of VCAM1 
expression. The microarray profile was successfully vali- 
dated by RT-PCR on six candidate genes (see Fig. 4) with 
the exception of Runx2 where an increase in Runx2 PCR 
product was detected at the TBL20 LPS 18 h time point. 
qRT-PCR validation indicated a 13.9-fold (log2) lower 
expression level for VCAMI in TBL20 versus BL20-LPS 
(Fig. S2, panels 4A), compared with a difference of 7.4- 
fold (log2) generated by the array data; while for Runx2, 
qRT-PCR showed a 5.1 -fold (log2) reduction in expression 
values between TBL20 and BL20-LPS (Fig. S2, panels 
4B), compared with a difference of 2.5-fold (log2) from the 
array data. Pathway analysis of this profile identified 
enrichment for genes within the categories 'inflammatory 
response' and 'inflammatory disease', and it was of inter- 
est that 'activation of IRF by cytosolic pattern recognition 
receptors' and 'interferon signalling' were identified as the 
most significantly enriched canonical pathways. 

Profile 8 is the reciprocal of profile 7 and represents 98 
genes where expression was downregulated by LPS 
treatment of BL20 but showed no change associated with 
infection of BL20 (BL20 versus TBL20). Like profile 7, 
genes sensitive to BW720c treatment were under- 
represented in this profile with only a single one present. 
The microarray data were validated by RT-PCR for three 
genes and by qRT-PCR for BCOR and RGS1 (Fig. 4 and 
Fig. S2, panels 8Aand 8B). Pathway analysis showed this 
profile is enriched for genes predicted to function in cell 
signalling, small molecule biochemistry and metabolism. 
This suggests that the infected cell maintains metabolic 
processes required for cellular proliferation (as also indi- 
cated for profile 6). 

Together, the eight profiles identified above provide evi- 
dence that the parasite-infected cell is associated with 
significant modulation of the response of uninfected BL20 
cells to activation with LPS. A sizeable proportion of these 
alterations to host cell gene expression are under control 
of a viable parasite, as they are altered by treatment with 
BW720C. However, the majority of the identified infected 
cell modulations to gene expression did not respond to 
drug treatment. While this may be influenced by the 
shorter time periods used for drug treatment relative to 
previous studies (Oura etal., 2006; von Schubert etal., 
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Fig. 6. Enrichment analysis of the infected cell-modulated gene list. The infected cell-modulated gene list (n= 1214) was analysed for 
evidence of enrichment with respect to two areas within the pathway analysis database, i.e. molecular function and canonical pathways. 
The gene list was also analysed for enrichment with respect to a number of denoted custom pathways. Top-ranking categories are presented 
together with the significance of the enrichment [-log(P-value)]. 



2010), it seems likely that a large number of the altera- 
tions to host cell gene expression are genuinely non- 
responsive. Such alterations could be derived from 
epigenetic events that occurred during establishment or 
passage of the BL20 and TBL20 cell lines. These could 
be disparate events that are independent of infection, 
although it is also possible that epigenetic events are 
influenced by the parasite to generate infected cell- 
associated alterations that cannot be altered by drug 
treatment over 48 h. 

Pathway analysis of the infected cell-modulated 
data set can predict phenotypes displayed by 
T. annulata-zhfecfec/ leucocytes 

Derivation of the eight profiles of genes displaying altered 
expression levels in infected TBL20 cells relative to LPS- 
stimulated uninfected BL20, provides evidence that 
T. annu/afa-infected cells alter the gene expression 
outcome associated with an activation event, i.e. the 
response to LPS. To gain information on the potential 
impact of these alterations, the entire list of 1214 genes 
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was analysed for enrichment in specific molecular func- 
tions and canonical pathways. A large number of catego- 
ries were identified for the 976 genes mapped onto the 
pathway analysis framework but it was clear that the data 
set was very significantly enriched for genes predicted to 
function in the inflammatory response, cell death, cancer 
and the antimicrobial response (Fig. 6). 

To make the output of the pathway analysis more pre- 
dictive, further enrichment analysis was undertaken, 
noting whether each molecule promoted (pro) or inhibited 
(anti) a particular molecular function. The proportion of 
genes designated as 'pro' or 'anti' for each functional 
category was then examined in relation to whether their 
expression is repressed or enhanced by infected TBL20, in 
order to infer how the infection-associated modulation 
could influence a functional outcome. Interestingly, func- 
tional categories highlighted by this analysis (Table 1 ) pre- 
dominantly show an association with a phenotype that is 
known or could be reasonably predicted for the infected 
cell. For example, in the 'metastasis of cells' subcategory, 
the observed trend indicated repression of anti-metastatic 
genes and elevated expression of pro-metastatic genes. 
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Table 1. 


Functional classification of the infected cell-modulated gene set. 












Functional classification 


Repressed in parasite-infected cell 


Enhanced in parasite-infected cell 


No. genes (pro/anti) 


P-value 


No. genes (pro/anti) 


P-value 


■| 


Inflammatory response 
















Activation of leucocytes 


46 


(27/9) 


2.7 X 1 0-'' 


13 


(9/1) 


3.8 X 10-3 




Antimicrobial response 


19 


(3/1) 


4.9 X 10"^ 










Activation of lymphocytes 


31 


(19/7) 


y.b X 1 U 




(5/1) 


1.2x10 




Antiviral response 


13 


(3/1) 


y.b X lU 


2 


(0/0) 


1.2x1 0-^ 




Inflammation 


30 


(1 2/3) 


/ X 1 u 


12 


(3/2) 


D.O X 1 u 


o 

d. 


Cell death 
















Apoptosis of tumour cell line 


72 


(34/22) 


2.3 X 10-^ 


20 


(12/7) 


3.1 X 10-3 




Apoptosis 


134 


(70/53) 


2.1 X 10"^ 


42 


(24/1 9) 


1.1 X 10-* 




Cell death of tumour cell line 


78 


(38/28) 


^.1 X 1 U 


24 


(13/9) 


y.4 X 1 u 




Cell death 


151 


(75/62) 


b.b X 1 u 


49 


(26/22) 


O O i A-5 
O.lL X 1 U 


3 


Cellular movement 
















Infiltration of cells 


29 


(10/6) 


4.2x10-^ 


- 




- 




Invasion of eukaryotic cells 


- 




- 


14 


(10/3) 


5.9 X 10-* 




Migration of antigen presenting cells 


14 


(7/1) 


1.6 X 10-"^ 






- 




Influx of leucocytes 


7 


(4/1) 


2.0 X 10"^ 


— 








Chemoattraction of leucocytes 


6 


(6/0) 


^.c. X 1 L) 








A 


Cancer 
















Metastasis of cells 


10 


(2/6) 


3.2x10-= 


13 


(3/1) 


2.3 X 10-3 




Tumorigenesis 


179 


(26/25) 


1.7 X 10-3 


60 


(6/2) 


1.7 X 10-= 




Cancer 


173 


(16/13) 


7.6 X 10-= 


56 


(3/1) 


1.6 X 10-= 




Disease of tumour 


18 


(2/4) 


7.3 X 10-* 










Growth of tumour 


20 


(2/9) 


1.2x10-3 










Tumour 








43 


(4/2) 


1.1 X 10^ 


C 

o 


Cell assembly and organization 
















Formation of actin filaments 








8 


(3/1) 


1.3x10-3 




Formation of actin cytoskeleton 


4 


(0/0) 


1.3 X 10-3 


4 


(2/0) 


1.1 X 10-^ 




Formation of filaments 








9 


(5/1) 


4.1 X 10-3 




Formation of cellular protrusions 








10 


(4/0) 


2.7 X 10-3 



implying that the parasite-infected leucocyte promotes 
metastasis. Genes classified within 'cancer' and 'tumori- 
genesis' were associated with a bias towards pro function 
in genes showing enhanced expression associated with 
the infected cell, whereas a bias towards anti-function was 
associated with repression by infected TBL20 within the 
categories 'disease of tumour' and 'growth of tumour'. 
Additional categories include the inflammatory response 
and cell movement. In the inflammatory response cat- 
egory, the data indicated that genes promoting an 
'anti-microbial/anti-viral' response or 'inflammation' are 
associated with repression by parasite-infected cells, while 
for 'activation of leucocytes' and 'activation of lymphocytes' 
a trend towards elevated expression of positive effector 
genes was obtained. Similar results were observed for cell 
movement, one subcategory 'invasion of eukaryotic cells' 
indicating elevated expression of genes pro this function, 
whereas the next three subcategories all indicated a trend 
towards repression of genes pro movement, including 
those categorized in 'influx of leucocytes' and 'chemoat- 
traction of leucocytes'. One possible explanation for these 
results is that they reflect a requirement for activation and 
dissemination of the infected leucocyte together with 
dysregulation/evasion of the immune response. 



The ability of the pathway analysis to predict pheno- 
types displayed by the 7he//er/a-infected cell suggests 
that the infected cell-modulated gene list could be mined 
to highlight candidate novel pathways and genes that 
operate in its establishment and maintenance. Examples 
of several candidate genes and categories of interest are 
given below. 

Sixty genes (24%) were placed in the cancer category 
(P-value, 1.56x10"^), including pro-cancer genes 
MMP-iS, SPARC, NDRG1, PAK1, ICAM1, ANGPT4, 
IL2RB-\\ke, RGS1, SELL, SERPINE2, CCR7, CTTN, 
LTB4R, SLC2A1 and 7P73 which are expressed at sig- 
nificantly higher levels in parasite-infected TBL20 relative 
to BL20-LPS. Intercellular adhesion molecule-1 (ICAM1) 
is elevated sixfold in TBL20 cells relative to BL20-LPS 
(profile 1), and elevated 23-fold relative to uninfected 
BL20. ICAM1 is a trans-membrane glycoprotein and is 
encoded by a classical NF-kB target gene (Huang etal., 
2004; Lin etal., 2006). It plays an important role in the 
immune response, but has also been implicated in cancer 
metastasis and linked to tumour prognosis and progres- 
sion (Huang etal., 2004; Lin etal., 2006). Interestingly, 
the level of ICAM1 mRNA was found to be significantly 
higher in infected monocytes derived from disease- 
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susceptible Holstein cattle relative to resistant Saliiwal 
(Jensen etal., 2008). Such a difference could be gener- 
ated by variability in the level of activated NF-kB in 
infected cells from different breeds, or a difference in a 
parasite-dependent mechanism that controls NF-kB 
target gene expression. Further investigations are 
required to characterize the role of ICAM1 expression in 
determining the virulence of the T. annu/afa-infected cell. 

Secreted protein acidic and rich in cysteine (SPARC) is 
a matricellular protein that influences cell-matrix interac- 
tions and has been reported to have multiple roles in 
macrophage function, inflammatory processes, cell mor- 
phology, inhibition of cell cycle progression and synthesis 
of the extracellular matrix (Rempel etal., 2001). IVIany 
cancer types exhibit increased SPARC levels upon inva- 
sion or metastasis (Rempel etal., 2001; Podhajcer etal., 
2008) and upregulation of SPARC expression (53-fold) in 
infected in TBL20 cells (profile 4) may play an important 
role in establishing the infected cell phenotype. For 
example, recent studies have described a connection 
between SPARC and the p53 tumour suppressor 
(Fenouille etal., 2011). Essentially, SPARC acts as an 
anti-stress factor by mediating degradation of p53 through 
AKT-mediated MDM2 phosphorylation to promote 
melanoma cell survival. Array data indicate that the level 
of p53 mRNA is high in 7A7e/7eria-infected leucocytes, but 
protein levels are virtually undetectable (B. Shiels and S. 
McKellar, unpublished). The PI3K/AKT pathway is known 
to be constitutively activated in the 7/?e//er;a-infected leu- 
cocyte (Baumgartner etal., 2000; Heussler etal., 2001). 
Investigation of a mechanism that may operate to 
promote constitutive degradation of p53 is required. 

The category 'cellular assembly and organization' was 
predicted to show an association with the infected cell 
phenotype, as a bias towards enhanced expression of 
pro genes within a number of subcategories was identi- 
fied. Subcategories included 'formation of actin filaments' 
and 'formation of actin cytoskeleton' and predict that 
infection of the host cell is associated with manipulation 
of cell shape and the actin cytoskeleton. Upregulated 
genes included CTTN, SDC4, ICAM1, PAK1, ELN, 
SPARC, SORBS1, TNS1 and F2R. Theileria infection 
is known to alter the morphology of the host cell 
(Schmuckli-Maurer etal., 2010; Baumgartner, 2011) and 
it has been postulated that organization of the actin 
cytoskeleton has an impact on parasite-dependent acti- 
vation of NF-kB via the IKK signalosome (Schmuckli- 
Maurer etal., 2010). Recent studies have shown that 
parasite-dependent polarization of actin-based structures 
such as lamellipodia and podosomes are critical for both 
attachment to substrates and motility (Baumgartner, 
2011). Molecules thought to function in regulation of 
these structures include the Src kinase, Hck, the Rho 
kinase, ROCK and cortactin. Elevated expression of the 
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cortactin (CTTN) and Hck (haemopoietic cell kinase) 
genes in infected TBL20 cells support a role in regulation 
of infected cell shape. 

The gene, PAK1, encoding the p21 (Cdc42/Rac)- 
activated kinase 1 (member of a family of serine/threonine 
protein kinases) is significantly upregulated (> 22-fold) in 
777e//er/a-infected TBL20 compared with BL20-LPS 
(profile 4). PAK1 has been reported to have a number of 
functions (e.g. cell survival, motility, proliferation) and is 
implicated in transformation and tumour progression 
(Huynh etal., 2010). PAK1 is a key molecule in the regu- 
lation of the actin and microtubule cytoskeleton (reviewed 
by Bokoch, 2003), plays a crucial role in formation of 
podosomes and membrane ruffles, can phosphorylate 
cortactin to regulate the dynamics of branched actin fila- 
ments (Webb etal., 2006) and enhances cell migration 
and proliferation viaAKT (Huynh etal., 2010). Intriguingly, 
PAK1 has been reported to operate in pathogen- 
dependent activation of NF-kB via NIK (Neumann etal., 
2006) and inhibition of NIK reduces activation of NF-kB in 
777e/7er/a-infected cells by 25% (Heussler etal., 2002). 
Further studies investigating the role PAK1 may play in 
establishing the transformed phenotype of the Thellerla- 
infected cell are indicated. 

In contrast to the categories and genes highlighted 
above, an expected bias towards upregulation of anti- 
apoptotic genes together with downregulation of pro- 
apoptotic genes could not be identified, although the 
present analysis framework may not have had the power 
to resolve such a phenomenon. This may be due to 
unequal numbers of 'pro' and 'anti' apoptosis genes 
denoted in the database together with the fact that 
genes functioning in apoptosis commonly switch from a 
'pro' to 'anti' function depending on cellular context (Fan- 
njiang etal., 2003; Hess etal., 2004) and such subtle- 
ties will not be reflected in the database. Moreover, the 
probability of progressing towards apoptosis may be 
associated with qualitative differences, such as phos- 
phorylation status, or the stoichiometry of a number of 
key factors and this would not have been revealed by 
the analysis. 

Theileria-znfecfed TBL20 cells modulate expression 
outcomes associated with activation of NF-kB and 
other Infection-activated transcription factors 

Stimulation of cells with inflammatory mediators such as 
LPS, and infection by T. annulata, is known to result in 
activation of transcription factors via a number of signal- 
ling pathways. Therefore, the infected cell-modulated data 
set of 1214 genes was analysed with respect to canonical 
pathways, including those involved in signalling. Unsur- 
prisingly, the results highlighted modulation of signalling 
events associated with the inflammatory response. Sig- 
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Fig. 7. NF-kB signalling pathway. The genes which are coloured are differentially expressed in TBL20 compared with BL20; red = up in TBL, 
green = down in TBL. The arrows are beside genes present in the infected cell-modulated list (1214). For example, -V^ '^ means super-up in 
BL20-LPS I up in TBL20 compared with BL20. The data for IL-1R/TLR and growth factor receptor are exploded in two boxes. 



nificant pathways include TNFR2, pattern recognition 
receptors, CD40, NF-kB, PI3K, LPS-stimulated IVIAPK 
and interferon signalling (Fig. 6). An event connmon 
among many of these pathways is the activation of NF-kB. 
As predicted from the findings of Heussler etal. (2002), 
the microarray data generated little evidence for major 
repression of the NF-kB signalling pathway below the 
level of the activated signalosome (Fig. 7); repressed 
expression of the genes encoding three IkBs (relative to 



LPS) would promote NF-kB activation but there is no 
evidence of a difference in expression of these IkB genes 
between non-stimulated BL20 and TBL20 cells. Upstream 
of the signalosome, repression of genes encoding recep- 
tors appears to be in operation and it can be concluded 
that, in general, infected TBL20 down regulate the ability to 
recognize extrinsic stimuli that induce activation of NF-kB. 
Exceptions to this observation are upregulation of the 
genes encoding the TNFR2 and TGFBR3 receptors. Both 
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TNF and TGF|3 signalling pathways have been associated 
with the 7/7e//er/a-infected leucocyte (Guergnon etai., 
2003; Ghaussepied etal., 2010) and could potentially 
activate NF-kB viaTRAF2 or PI3K. However, evidence for 
a role of the PIK3 pathway in activating NF-kB is conflict- 
ing (Baumgartner etai., 2000; Heussler etal., 2001), and 
the importance of a TNF autocrine loop is unclear (Guer- 
gnon etal., 2003). Furthermore, microarray data indicate 
that the gene encoding tumour necrosis factor alpha 
induced protein 3 (TNFAIP3/A20) that acts as an inhibitor 
of TNFR2-dependent NF-kB signalling (Charan etai., 
2011) is upregulated in TBL20 relative to BL20. The 
importance of the PI3K pathway for proliferation of the 
infected leucocyte (Baumgartner etal., 2000; Heussler 
etal., 2001) is supported in the present study by evidence 
for upregulation of the genes encoding PI3K regulatory 
subunits 5 and 6 relative to repression by the LPS 
response. 

Using a series of custom pathways, the infection- 
modulated data set was also screened for enrichment of 
target genes of transcription factors activated by the 
parasite-infected cell and associated with the inflamma- 
tory response. As shown in Fig. 6, significant enrichment 
of NF-kB targets and cMyc repression targets was 
observed, together with a particular enrichment of genes 
regulated in response to stimulation by TGF(3. This pro- 
vides evidence that the T/?e//eria-infected leucocyte sig- 
nificantly alters the outcome of cellular activation by 
modulating the expression of target genes regulated by 
constitutively activated transcription factors. 

Investigation of NF-kB target genes revealed that the 
majority (46) were repressed in infected TBL20 cells rela- 
tive to LPS stimulated uninfected cells (BL20-LPS), while 
15 showed evidence of enhanced expression (Fig. 8). For 
repressed NF-kB gene targets, enrichment was observed 
for pro-inflammatory response (e.g. CXCL10, CCL4 
and IL15), immunoreceptors and pro-differentiation of 
leucocytes. CXCL10 is an angiostatic chemokine that 
represses the ability of Burkitt's lymphoma cells to 
produce tumours in nude mice (Sgadari etai., 1996). 
Overexpression of CXCL10 has been implicated in inhi- 
bition of the proliferation process in a variety of cell types. 
IL15 is a major pro-inflammatory cytokine involved in 
enhancing cytotoxic CDS T cell and NK cell activity and 
has therapeutic indications against cancer (reviewed by 
Steel etai., 2011). Based on the proposed role of CDS T 
cells and NK cells in the immune response against T. an- 
nulata (reviewed in Preston etal., 1999), major down- 
regulation of IL15 by the parasite-infected cell (i.e. 100- 
fold less in TBL20 compared with BL20 LPS) could be 
predicted to be beneficial to its establishment. Repressed 
NF-kB target genes also included regulators of apoptosis 
and genes involved in regulating the NF-kB mediated 
inflammatory response, i.e. TRAF1 (Henkler etal., 2003), 
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CFMR/FLIP (Zaitseva etai., 2011) and CASP4 (Laksh- 
manan and Porter, 2007; IVIao etai., 2010). 

The number of infected cell-associated upregulated 
NF-kB target genes (Fig. S) was too small for enrichment 
analysis, but did include some interesting members, such 
as genes encoding the pro-metastatic chemokine recep- 
tor CCR7 (Muller etal., 2001), ICAM1, SPARC (see 
above) and the N-myc downstream regulated 1 (NDRG1). 

Investigation of genes classified as repression targets 
of cMyc indicated that the majority, 7 out of 10, displayed 
an expression profile consistent with maintenance of 
repression in the infected cell (data not shown). This 
included two genes encoding inhibitors of cell cycle pro- 
gression, CDKN1A and CDKN2B (Cartel and Shchors, 
2003). In marked contrast, NDRG1 (classified as regu- 
lated by NF-kB and cMYC repression) showed a signifi- 
cant elevation associated with the parasite-infected cell 
relative to repression in BL20-LPS (Fig. 8). NDRG1 has 
been proposed to be a suppressor of tumour growth and 
an attenuator of NF-kB activation (Hosoi et ai., 2009) and, 
therefore, elevation of expression by infection could be 
considered counter-intuitive. However, a role for NDRG1 
at the spindle checkpoint of mitosis has been established, 
and NDRG1 is thought to be important for maintaining 
correct microtubule functioning (Kim etal., 2004) and 
spindle formation. Recent work has shown an intricate 
relationship between the Theileria schizont and host cell 
mitotic and central spindles, dependent on recruitment of 
host cell Polo-like kinase (von Schubert etai., 2010). 

Evidence for modulation of TCPp target genes was also 
obtained, although a uniform alteration of the predicted 
response to TGFp stimulation could not be identified for 
TBL20 cell relative to BL20-LPS (data not shown). Thus, 
roughly equivalent numbers of target genes expected to 
be upregulated or downregulated by TCPp were identified 
within the 'infected cell-associated repressed' gene list, 
with a similar finding for the 'enhanced' gene list. The data 
set of TGFp target genes modulated by infected TBL20 
was found to be enriched for genes placed in the catego- 
ries of 'cancer', 'cell death', 'cellular growth' and 'cellular 
development'. Of particular note was that TGFp targets 
categorized as inflammatory response genes were 
generally associated with repression of expression by 
parasite-infected TBL20 relative to the BL20-LPS 
response. The lack of any association between a pre- 
dicted TGFp response and alteration of expression by the 
infected cell may be explained, in part, by the finding that 
expression of genes encoding TGFP2 and TGFps and the 
TGFpRS receptor is slightly elevated in parasite-infected 
TBL20 cells, while TCPpi is slightly reduced. The TGFp 
signalling pathway has been implicated as a determinant 
of pathogenicity and metastatic activity in tropical theile- 
riosis: reduced TGFp2 expression was measured for 
infected cells derived from an attenuated vaccine line 
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Fig. 8. Hierarchical clustering of NF-kB response genes. The NF-kB targets genes identified in the infected cell-modulated gene list were 
selected for hierarchical clustenng analysis and the results are presented as a heat map. Genes highlighted in bold show altered array 
expression following BW720c treatment for 48 h to kill the parasite. Underlined genes show altered RT-PCR expression following BW720c 
treatment for 72 h. 
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and disease-resistant Sahiwal (60s indicus) cattle 
(Ciiaussepied etal., 2010). Differential expression of 
TGFp target genes was also identified, with elevated as 
well as repressed expression reported for the attenuated 
vaccine cell line. Thus, modulation of TGFp signalling and 
the respective target genes appears to be a common 
event associated with T. annulata-'miected leucocytes. 

The results of this study provide evidence that the 
T. annu/afa-infected cell tightly controls the outcome of 
NF-kB activation. This modulation, in general, can be 
predicted to be beneficial to establishment of the infected 
cell. The data also indicates that the outcome of activating 
additional pathways/transcription factors l<nown to be 
associated with TtoVer/a-transformed leucocytes are 
modulated, e.g. TNF, TGFp and cMyc. Thus, a major 
re-organization of networl<s that regulate gene expression 
changes associated with the inflammatory response and 
cancer is lil<ely to be concomitant with transformation of 
the infected leucocyte. 

Recent studies on the related apicomplexan parasite, 
T. gondii have identified parasite proteins that regulate 
inflammatory gene networks by activating transcription 
factors, including NF-kB (Saeij etal., 2007; Rosowsl<i 
etal., 2011). Evidence has also been generated that 
parasite-dependent modification of histone H3 can 
prevent transcription factors binding to host-cell inflamma- 
tory gene promoters (Long and Denkers, 2009). It has 
been proposed that such a mechanism could account for 
the large-scale downregulation of pro-inflammatory genes 
by T. gondii (Leng etal., 2009). While it is possible that 
Theileria parasites deploy similar mechanisms to manipu- 
late gene expression on a wide scale, our results suggest 
additional alterations are involved. This is indicated by the 
four pairs of reciprocal modulated expression profiles 
(Figs 3 and 4), the ability of infected cells to reverse gene 
expression changes induced by LPS and the non- 
responsive and non-reversible nature of a large number of 
the host cell gene expression changes when the parasite 
is killed by drug. Such changes could be generated by 
an irreversible reconfiguration of transcription factor/ 
expression networks associated with differentiation 
events, as increasingly proposed for stem cell systems 
(Chickarmane etal., 2009; Wang etal., 2010). In addition, 
deployment by T. annulata of DNA-binding proteins to the 
host cell nucleus that show similarity to mammalian 
HIVIGA proteins could be involved (Swan etal., 1999). 
HMGAs act as hubs of nuclear function and play a major 
role in organization of gene expression networks associ- 
ated with inflammation and cancer (reviewed in Cleynen 
and Van de Ven, 2008). Given the propensity of Theileria 
parasites to control gene expression networks associated 
with inflammatory disease and cancer, further study of the 
mechanisms deployed and the outcome they promote is 
warranted. 
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Experimental procedures 

Cell culture and LPS stimulation 

Low-passage BL20, an uninfected bovine lymphosarcoma cell 
line (Olobo and Black, 1989), and TBL20, a T. annulata (Ankara 
strain) infected BL20 cell line (Shiels etal., 1986), were cultured 
in standard complete medium (Swan etal., 1999) except that 
heat-inactivated fetal bovine serum (FBS; Sigma) was used. 
Cultures were seeded with 2x10^ cells mM and maintained by 
feeding with fresh medium every 2-3 days (Schmuckli-Maurer 
etal., 2010). Stimulation of cells was carried out for up to 48 h in 
six-well plates with 5x10'^ cells welM in a total of 5 ml of 
medium, set up immediately prior to addition of 1 |ig mM LPS 
(Sigma: L2630). Optimum LPS concentration was determined by 
preliminary titration experiments which measured the minimum 
concentration of drug required to stimulate cells, based on NF-kB 
reporter activity, without causing rapid cytotoxicity (data not 
shown). LPS was dissolved in dimethyl-sulfoxide (DMSO) and 
control cells were incubated with DMSO alone. Treatment of 
TBL20 with BW720C was performed as described previously 
(Schmuckli-Maurer etal., 2010). 

Cell viability counts and apoptosis assay 

Assessment of viable and non-viable cell numbers was per- 
formed using Trypan Blue exclusion counted using a haemocy- 
tometer; viable cells exclude the dye. Three replicate cultures of 
BL20 and TBL20 were seeded with 2.5 x 10^ cells mr' and incu- 
bated at 37°C for 18 h in the presence of LPS (1 |xg mr') or 
DMSO (1 |il mf) as no drug control. Control experiments com- 
paring cells culture in the presence or absence of DMSO showed 
no evidence of an effect on cell viability or activation of NF-kB 
(data not shown). Caspase-3/7 activities were measured using a 
Caspase-3/7 Assay kit (Promega) according to the manufactur- 
er's protocol. Each assay was performed using 2.0 x 10"* cells 
from three replicate cultures per experimental condition at 18 h. 
Luminescence was measured using a TD -20/20 luminometer 
(Turner Designs). A no-cell control was carried out to account for 
baseline signal generated from fresh complete medium. 

NF-KB-dependent reporter assay (dual luciferase assay) 

Estimation of NF-kB activity was performed by transient 
co-transfection of a reporter plasmid containing a promoter with 
an NF-KB3 recognition sequence linked to Firefly luciferase (NF- 
kB-LUC) and a Renilla luciferase plasmid (phRG-TK renilla), 
followed by sequential measurement of Firefly and Renilla luci- 
ferase activity, as previously described (Schmuckli-Maurer etal., 
2010). Two micrograms of NF-kB3-Luc and 0.5 ng of phRG-TK 
Renilla plasmids were mixed with 4x10'' BL20 or TBL20 cells in 
300 ^il of RPMI-1 640-2% FBS without antibiotics and electropo- 
ration performed using a Gene Pulser (Bio-Rad) set at 220V, 
950 [iF. After electroporation, cells were immediately added to 
5 ml of pre-warmed complete medium in six-well plates. After 
24 h recovery, LPS (1 |xg mf) or DMSO were added for 4 h. 
Samples were harvested and the cell pellet washed with PBS 
(phosphate buffer saline). 

A Dual Luciferase® Reporter Assay System (Promega) was 
used to assay activity, using the manufacturer's protocol. NF-kB 
promoter activity was calculated by normalizing firefly luciferase 
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activity relative to the activity of constitutive Renilla luciferase. 
Each electroporation and luciferase assay was performed in 
quadruplicate, 20 ^1 of cell extract was used per assay. 

Statistical analysis 

Statistical analysis and generation of graphs were performed 
using GraphPad Prism Software version 5.0 (La Jolla, CA). Data 
were expressed as mean ± standard error of the mean (SEM). 
Differences and significance between mean values were calcu- 
lated using a two-tailed Student's f-test with P-values £ 0.05 
considered to be statistically significant. 

Indirect immunofluorescence analysis 

Indirect immunofluorescence analysis (IFA) was performed on 
uninfected BL20 and infected TBL20 cells cultured with LPS or in 
DMSO (1 \i\ mM), as no drug control, for 4 h. Cytospin slides 
were generated using poly L-lysine glass slides and IFA per- 
formed, as described previously (Schmuckli-Maurer ef a/., 2010). 
Primary antibodies, anti-IKKy/NEMO (BD Biosciences; 611306), 
anti-NF-KB p65 (Santa Cruz; sc-372) and anti-NF-KB p65 (Santa 
Cruz; sc-8008), were used at 1:50 dilution; anti-rabbit IgG and 
anti-mouse IgG secondary antibodies, conjugated to Alexa 488 
(Invitrogen; A-11034; A-11029), were used at 1:200 dilution. 
Images were acquired using an Olympus BX60 microscope, 
SPOT camera and SPOT™ Advanced image software Version 
Mac: 4.6.1.26 (Diagnostic Instruments). 

Western blotting 

Total cell extracts were prepared from approximately 5x10^ cells 
grown to a density of 1 x 10" cells mh' in a six-well plate, resus- 
pended in 100|xl of PBS and lysed with an equal volume of 
SDS sample buffer (2x concentrate, Sigma; S3401). Extracts, 
adjusted for equal protein loading, were resolved by SDS-PAGE 
(Laemmli, 1970) with subsequent immunoblotting and detection 
as described previously (Schmuckli-Maurer etal., 2010). The 
anti-TLR4 antibody (Santa Cruz-Biotechnology, sc-10741) was 
used at 1:1000 dilution and anti-tubulin antibody (Sigma, T9028) 
at 1 :5000. 

RNA isolation and processing 

Total RNA for analysis by microarray was isolated in triplicate 
from 30 ml of cultures of infected and uninfected cells and cells 
treated with LPS for 4 and 1 8 h. Cell pellets of 1 .5-2.0 x 1 0' cells 
were resuspended immediately in TRI Reagent (Sigma; T9424), 
and extracted according to the manufacturer's instructions. To 
eliminate possible contamination with genomic DNA, RNA 
samples used for RT-PCR and QRT-PCR were treated with 
DNase I (Qiagen; 79254) either before purification on a silica-gel 
membrane-based RNeasy® column (Qiagen; 74104) or by an 
on-column DNase I treatment. A2m/A2so ratios obtained for all 
samples fell in the acceptable range of 1.78-2.1 and quality and 
integrity of the total RNA preparations was confirmed by gel 
electrophoresis. 

Custom bovine microarray and analysis 

The oligonucleotide microarray used in this study included all 
bovine RNA RefSeq sequences available at NCBI (26 751) cor- 



responding to the Btau 4.0 assembly of the bovine genome 
(http://www.ncbi.nlm.nih.gov/RefSeq/) and a number of putative 
bovine RNA Sequences and ESTs (4730) kindly provided by 
Kirsty Jensen, Roslin Institute, University of Edinburgh. Only 
extant sequences in the RefSeq mRNA database (December 
2011) or UMD3.1 genome assembly were used in the analysis 
and this represents 19 777 sequences. Design of oligonucle- 
otides, array synthesis, hybridization to the array and normaliza- 
tion of array data were carried out by Roche NimbleGen, 
Madison, USA. Each gene was represented by two identical sets 
of five 60-mer oligonucleotide probes which were isothermal with 
respect to melting temperature. cDNA was generated from 1 0 |j,g 
of total RNA using oligo(dT) primer and tagged with 3'-Cy3 dye 
after which 13 |ig of labelled cDNA was hybridized to the array. 
Gene expression values were calculated from a RMA-normalized 
data set (Irizarry ef a/., 2003) and differentially expressed genes 
were identified using Rank Product Analysis (Breitling etal., 
2004). Genes were defined as differentially expressed using the 
criteria of a false discovery rate of less than 5% and a fold change 
of greater than two. In addition to the bovine sequences, the 
array contained oligonucleotide probes representing the T. annu- 
lata genome, although these data were not analysed in the 
course of the present study. The data discussed in this publica- 
tion have been deposited in NCBI's Gene Expression Omnibus 
(Edgar etal., 2002) and are accessible through GEO Series 
Accession No. GSE36428 (http://www.ncbi.nlm.nih.gov/geo/ 
query/acc.cgi?acc=GSE36428). 

To assess the relationship between the differentially expressed 
gene sets identified by RPA in TBL20 and LPS stimulated BL20 
cultures, a Monte Carlo simulation was used to estimate the 
overlap expected to occur by chance (Metropolis and Ulam, 
1949; Kabakchiev etal., 2010). A chi-squared test was used to 
statistically test whether the observed overlap was greater than 
that expected by chance. 

In order to identify gene expression networks that were 
enriched for genes highlighted in this study. Ingenuity Pathway 
Analysis V9 (Ingenuity® Systems, http://www.ingenuity.com) was 
utilized when possible. Each bovine RNA sequence represented 
on the array was mapped to its corresponding gene in the Inge- 
nuity Pathways Knowledge Base. The list of canonical pathways 
in the system was augmented with a number of specific path- 
ways, namely a list of 453 validated and predicted NF-kB target 
genes (Sharif etal., 2007), the NF-kB pathway the cMyc 
pathway, validated targets of cMyc transcriptional repression, 
validated targets of cMyc transcriptional activation, AP-1 target 
genes, ATF2 network genes and the WNT signalling pathway, 
largely derived from information in the NCBI BioSystems data- 
base. The significance of the association between novel gene 
lists and each pathway was measured in two ways. First, a ratio 
was calculated for the number of genes in a given list that map to 
the pathway divided by the total number of genes on the array 
that map to the same pathway. Second, Fischer's exact test was 
used to calculate a P-value to determine the probability that an 
association between a given gene list and a particular canonical 
pathway occurred by chance alone. 

Reverse transcription PCR (RT-PCR) 

Relative expression levels of selected candidate genes of inter- 
est were analysed by semi-quantitative PCR using Superscript 
Ill-One Step RT-PCR (Invitrogen, 12574-026) and the standard 
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manufacturer's protocol. Primers were designed and tested for 
specificity using NCBI's online program (http://www.ncbi.nlm. 
nih.gov/tools/primer-blast/index.cgi) and oligonucleotide pro- 
perties calculator (http://www.basic.northwestern.edu/biotools/ 
oligocalc.html). Primer sequences, GenBank accession num- 
bers, annealing temperatures and product lengths for all gene- 
specific primers used in this study are given in Table S9. Primers 
were synthesized by Eurofins MWG Operon (Ebersberg, 
Germany) and whenever possible one of the primer pair spanned 
an intron/exon junction so that only cDNA was amplified or the 
primer pair amplified across at least one intron to give differen- 
tially sized genomic and cDNA products. Initial PGR optimization 
reactions were performed to ascertain the best annealing tem- 
perature (Tm), magnesium concentration and overall PGR effi- 
ciency. RT-PGR was performed in a final reaction volume of 25 |il 
with gene-specific primers at 200 nM and 12.5 ng of RNA tem- 
plate. Thermal cycling conditions were: an initial cycle of 50°G for 
30 min, 94°G for 2 min, followed by 36 cycles of 94°G 15 s, X°G 
for 30 s (where X°G indicates the experimentally determined 
annealing temperature) and 72°G for 1 min; a final extension step 
was performed at 72°G for 10 min. 

RT-PGR products were resolved relative to a 1 00 bp DNA 
ladder by electrophoresis through ethidium bromide-stained 1% 
TAB agarose gels. For all semi-quantitative RT-PGR assays, the 
expression levels of candidate genes were normalized to the 
levels obtained with primers for the genes encoding Beta 
actin (ACTp) or glyceraldehyde-3-phosphate dehydrogenase 
(GAPDH). These constitutively expressed genes did not show 
significant changes on the present array experiment and have 
been used previously as normalizing genes for bovine gene 
expression studies in T. annulata-'miected cells (Sager etal., 
1997). 

SYBR QRT-PGR methodology was utilized to quantify more 
accurately the relative fold change in mRNA expression for 
selected differentially expressed genes. cDNA synthesis was 
carried out using 2 \ig of DNase l-treated total RNA and an 
Affinityscript qPGR cDNA synthesis kit (Stratagene, USA; 
600559) with oligo(dT) primers, according to the manufacturer's 
instructions. Gontrol reactions without reverse transcriptase were 
always included. Real-time florescence detection was performed 
using 1 |xl of cDNA with Brilliant SYBR Green qPGR Master Mix 
buffer (Stratagene, 600548) according to the manufacturer's pro- 
tocol. The final concentration of each primer was 200 nM. A 
Stratagene Mx3005P Real-Time PGR System was used with the 
following thermal cycling parameters: 10 min at 95°G; 40 cycles 
of 30 s at 95°G (denaturation), 30 s annealing at X°G (Tm of 
primer pairs minus 3°G) and 30 s at 72°G (extension). After 40 
cycles of amplification, a melting curve analysis was carried out 
to verify the correct product by its specific melting temperature 
(Tm). All qPGR data were captured and analysed by Stratagene 
MxPro v4.10 software. Quantitative values were normalized rela- 
tive to those obtained for housekeeping genes, either p-actin or 
GAPDH, and fold changes calculated relative to the calibrator 
(DMSO-treated uninfected BL20 cells) using -2-"^' equation 
(Livak and Schmittgen, 2001). Statistical analysis and generation 
of graphs were performed using a two-tailed Student's f-test with 
GraphPad Prism Software version 5.0 (La Jolla, GA). 
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Fig. S1. Immunofluorescence analysis of activated BL20 cells in 
response to LPS stimulation. Panels 1-3: IKK signalosomes 
were detected in BL20 and LPS stimulated BL20 using a mono- 
clonal anti-IKKy antibody. IKK signalosomes (green dots) were 
observed only in BL20 LPS. DAPI nuclear staining is shown in 



panel 3. The arrow denotes a cluster of IKK signalosomes. 
Panels 4-6: NF-kB was detected in BL20 and LPS stimulated 
BL20 using anti-p65 antibody. Elevated expression and nuclear 
localization of NF-kB was observed in BL20 LPS. DAPI nuclear 
staining is shown in panel 6. Bar= 15 \im. 
Fig. S2. Panel of qRT-PCR results. 

Fig. S3. Top ranking cellular and molecular functions from 
pathway analysis. 
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